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Abstract 

An algorithm for numerical computation of the barrier transparency for the poten- 
tials surrounding Schwarzschild black holes is described for massive scalar particles. 
It is then applied to calculate the total (including all particle species and the contri- 
butions of twisted field configurations connected with Dirac monopoles) luminosity 
for the pion-kaon Hawking radiation from a Schwarzschild black hole with mass 
M = 10 12 g. It is found that the contribution due to monopoles can be of order 
11 % of the total pion-kaon luminosity. 
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1 Introductory remarks 



In a number of papers [1,2] it was shown that owing to nontrivial I 2 x § 2 - 
topology of black hole spacetimes the topologically inequivalent configurations 
(TICs) of complex scalar field can exist near black holes. Physically the ap- 
pearance of TICs for the given field should be obliged to the natural presence 
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of Dirac monopoles on black hole and because of the interaction with them 
complex scalars splits into TICs though the total (internal) magnetic charge 
of black hole remains equal to zero (see for more details Refs. [1,2] and be- 
low). As a consequence, description of the Hawking radiation process should 
be modified to include the contributions from different TICs. 

The expression for the luminosity regarding Hawking radiation for any TIC, 
however, contains an element of the S'-matrix connected with some potential 
barrier surrounding black hole and which effectively arises for the quantum 
scalar particle leaving black hole. One should therefore solve some scattering 
problem on the whole axis with the given potential. The relevant potentials, 
however, depend on both the type of black hole and the particle masses and are 
complicated since they can be described only in an implicit form. Therefore for 
physical results to be obtained one needs to apply the numerical methods. In 
the massless case the scattering problems discussed are simpler to treat and 
we obtained numerical results for luminosities in the case of Schwarzschild 
(SW) and Reissner- Nordstrom (RN) black holes in Refs. [3]. The massive 
case, however, is more realistic but the corresponding scattering problems are 
more complicated. Recently they have been studied in Ref. [4] in the SW case 
where it was possible to find an algorithm for numerical calculation of the 
conforming S'-matrix element and, accordingly, of the conforming Hawking 
radiation luminosity for massive scalar particles. The present paper contains 
a description of the algorithm and applies it to calculate the total luminosity 
for pions and kaons (both neutral and charged ones) for a SW black hole 
with mass M = 10 12 g. Of course, the given scalar mesons are not genuine 
elementary particles since they are composed of quarks but one can assume 
that the quark hadronization to some extent occurs inside black holes. Actually 
such an assumption is reasonable in light of a recent solution to the puzzle of 
the spin content of the proton, which involves assuming the proton contains 
virtual pions (as well as virtual quarks, antiquarks and gluons) and then the 
observed asymmetry can arise quite naturally (see, for instance, Ref. [5] for 
more details and further references). It is therefore reasonable to assume that 
the more massive objects, such as black holes contain virtual light mesons, 
such as pions and kaons. As a result, the task under consideration here makes 
sense. 

We write down the black hole metric under discussion (using the ordinary set 
of local coordinates t, r, ip) in the form 

ds 2 = adt 2 - a- x dr 2 - r 2 (d$ 2 + sin 2 tidip 2 ) (1) 

with a = 1 — r g /r, r g = 2M and M is the black hole mass. 

Throughout the paper we employ the system of units with h = c = G = 1, 
unless explicitly stated otherwise. Also we shall denote the set of the modulo 
square integrable complex functions on any manifold F furnished with an 
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integration measure as L 2 (F). When computing we use the following values 
of the scalar meson masses (see Ref. [6]) 

fi (TT°) = 134.9764 Mev,/i (7T ± ) = 139.56995 Mev, 
fi (K°, K°) = 497.672 Mev, /i (^ ± ) = 493.672 Mev . (2) 

2 Description of algorithm 

As was disscussed in [1,2] TICs of a complex scalar field on black holes are 
conditioned by the availability of a countable number of complex line bundles 
over the M 2 x § 2 -topology underlying the 4D black hole physics. Each TIC 
corresponds to sections of a complex line bundle E, which can be characterized 
by its Chern number n G Z (the set of integers). TIC with n = can be 
called untwisted, while the rest of the TICs with n^O should be reffered to as 
twisted. Using the fact that all the mentioned line bundles can be trivilized over 
the chart of local coordinates (t, r, </?) covering almost the whole manifold 
R 2 x S 2 one can obtain a suitable wave equation on the given chart for TIC <fi 
with mass /i and Chern number n G Z in the form 

(\g\)- 1/2 (d, - ieA^dg]) 1 ^ - ieA v )4>\ = , (3) 

where vector-potential for the corresponding Dirac monopole is A = A^dx^ = 
— j cos fidip. As was shown in Refs. [1], integrating F = dA over the surface t = 
const, r = const with topology § 2 gives rise to the Dirac charge quantization 
condition 

/n 
F = 4ir- = Airq 

s 2 6 

with magnetic charge q, so we can identify the coupling constant e with electric 
charge. Besides, the Maxwell equations dF = 0, d*F = are fulfilled [1] with 
the exterior differential d = dtdt + d r dr + d^dd + d^dip in coordinates t, r, </?, 
where * denotes the Hodge dual form. Also it should be emphasized that the 
total (internal) magnetic charge Q m of black hole which should be considered 
as the one summed up over all the monopoles remains equal to zero because 

Qm = - Y] n = . 

As was shown in [1,2] Eq.(3) has in L 2 (M 2 x S 2 ) a complete set of solutions of 
the form 

fwmin = -e lut Y n i m (§, ^Ruinir), Z = jn.|, \n\ + 1, |m| < I, (4) 
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where the explicit form the monopole (spherical) harmonics 5^u m (i?, ip) can be 
found in [1]. 



We now consider the functions -R^„(r). Denoting k = r g u, y(x) = r/r g and 
introducing the functions ip(x, k, I, n) = Ru(k)in( r ( x )), where y(x) is a function 
reverse to x(y) = y + ln(y — 1) (i.e., — oo < x < oo, 1 < y < oo), we find that 
ij)(x,k,l,n) obeys the Schrodinger-like equation 



( al 2 \ 

[ -j^p + k 2 I ip(x, k, I, n) = V(x, I, n, fJ,)i/)(x, k, I, n) , 



(5) 



V{x,l,n,n) = 



y(x) 



2 /(/ + l)-n 2 
^ + - + 



1 



(6) 



y 2 (x) y 3 (x) 

where \i = /io r g, or > through usual units, if /i and M in g then \i = 2/iqM/ii 2 ^ 
while if fio in Mev and M in g then 



H = 2- 3.762426569 • 10 _1 VoM , 

while = ch/G is the Planck mass square in g 2 . It is clear that as x — > +oo, 
V(x, I, n, //) — > fi 2 and then condition, which holds in the massless case 



/ 



Vldx < oo 



is violated and we cannot pose the scattering problem for Eq. (5) as in the 
massless case [3]. As an example, in Fig. 1 the numerically computed potential 
V is shown for charged kaons at n = 2, M = 10 15 g. 

Fig. 1. Typical potential barrier for charged kaons. 



As was shown in Ref. [4], the correct statement of the scattering problem for 
Eq. (5) with potential (6) consists in searching for two solutions ip + (x, k + , I, n), 
i/)~(x,k,l,n) of the equation (5) satisfying the following conditions 



i+i £,+ 7 , _ f e ikx + s l2 (k,l,n)e ikx + o(l), x -> -oo, 

if> [x,k ,l,n) - < Sn ( fc)Z)n ) w i(_2iA;+x) +o(l), x -> +oo, 



-ifca; 



j S 2 2 (A;,/,n)e ^ + o(l), x - -oo, 

^ (£, ft, t, nj — j i (2i/c + x) + s 2 i(k, I, n)w_ ifl i(—2ik + x) + o(l), a; -> +oo, 

2 ' " (7) 

2 2 

where fc + (A;) = -Jk 2 — ji 2 , p, = jj^yi = fc2 ^_ M 2 and the functions w ±i ^i (±2) are 
related to the Whittaker functions W ±i ^i(±z) (concerning the latter ones see 
e. g. Ref. [7]) by the relation 



w 



±i^±z) = W ±l ^(±z)e-^ 2 
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so that one can easily gain asymptotics (using the corresponding ones for 
Whittaker functions [7]) 

w iflt i(-2ik + x) = e ** + v*H2* + *l[i + 0{\k + x\-% x^+oo, 

W^i{2ik + X) = e -ik+x e -mn\2k+x\^ + x ^ +OQ , (g) 

We can see that there arises some S- matrix with elements Sij,i,j = 1,2 and 
it satisfies certain unitarity relations (for more details see Ref. [4]). As can be 
seen, the equation (3) corresponds to the lagrangian 

where the overbar signifies complex conjugation, T>^ = — ieA^ and, as a 
result, we have the conforming energy- momentum tensor for TIC with the 
Chern number n 



TV = Re[(^0)(P^) - -g^iV^iV^) - . 

Having obtained these relations, one can discuss the Hawking radiation process 
for any TIC of complex scalar field. Actually, one should use the above energy- 
momentum tensor to find, according to the standard prescription (see Refs. 
[1,2] and cited therein) and employing the asymptotics (7)-(8), the luminosity 
L(n) of Hawking radiation for TIC of a massive complex scalar field with Chern 
number n and mass /i is (in ordinary units) 

L{n) = hm 4vrr 2 < T tr >= A £ (2/ + l) Q (n), cj(n) = J ^f^ (9) 

with the vacuum expectation value < T tr >, where A = ^ (^fj) ~ 0.273673- 
10 50 erg-s _1 • M~ 2 (M in g), while the barrier transparency T(k,l,n) = 
\s u {k,l,n)\ 2 . 

It is clear that for the all configurations luminosity L of a black hole with 
respect to the Hawking radiation concerning the complex scalar field with a 
mass fj, to be obtained one should sum up over all n 

L = £ L(n) = L(0) + 2 £ L(n) , (10) 

nGZ n=l 

since L(—n) = L(n). Obviously, for neutral scalar field there exists only L(0) 
since it does not interact with monopoles. Also it is evident that for numerical 
computation of all the above luminosities one needs to have some algorithm 
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for calculating Su(k,l,n). This can be extracted from the results of Ref. [4]. 
Namely 

s u (k, I, n) = 2ik/[f-(x, k, /, n), f + (x, k + , Z, n)\ , (11) 

where [,] signifies the Wronskian of functions f~,f + , the so-called Jost type 
solutions of Eq. 5. In their turn, these functions and their derivatives obey the 
following integral equations 



f-(x,k,l,n) = e~ rkx + - / sm[k(x-t)]V(t,l,n,/j,)f-(t,k,l,n)dt, (12) 



(f-)' x (x,k,l,n) = -ike~ ikx + J cos[k(x - t)]V(t, I, n, fJ,)f~(t, k, I, n)dt , (13) 

— oo 

f + (x, k + ,l,n) =w ijly i(-2ik + x) + 

+oo 

— J ImK A) i(-2 ? A; + x) W _^i(2 ? A; + t)]V + (t,/,n,/i)/ + (t,A; + ,/,n)rft, (14) 



(f + y x (x,k + ,l,n) = ^w^(-2ik + x)+ 



+oo 



J \m[j-w l ^_(-2ik + x)w_ l ^ik + t)]V + ^ (15) 

X 

where the potential 



V + (x, I, n, ji) = /i 2 



1 



x y(x) 



+ 



1 - 



y(x) 



Z(/ + l)-n 2 ^ 1 



y 2 (x) 



y 3 (x) 



tends to when x — > +oo. Eqs. (12)-(13) are suitable for direct calculation 
while this is not the case for Eqs. (14)— (15) — it is difficult to evaluate the 
functions w(z) in (14)-(15) since there is no effective method to compute the 
Whittaker functions W(z) at z required. In the present paper we, therefore, 
stick to the following strategy. We employ the asymptotic expressions of w- 
functions at x — > +oo from Eq. (8) and replace Eqs. (14)-(15) by 



f + (x, k + ,l,n) = e 



i(k+ x+fi\n{2k+ x)) 



+ 



— [ lm[e i ^ k+ ^- t ^ ln ^}V + (tJ,n,ii)f + (t,k + J,n)dt, 

X 

(f + )' x (x, k + ,l,n) = ie i{k+x+ ~^ 2k+x) \k + + fi/x) + 

/ ImWl + ^y(k + (*-t)+iiHx/t))} V +U : z> n> fAf+(t, k + , I, n)dt . 
J k + x 



(16) 



(17) 
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The latter equations are appropriate for numerical evaluation but under this 
approximation, as numerical experiment shows, reliable calculations are only 
possible for /i < 0.01. This upper limit corresponds to /i (7r ± ) and M ~ 10 13 
g, so that for kaons the calculations are unreliable for M = 10 13 g. Therefore 
in the present paper we restrict ourselves to considering pions and kaons for 
M = 10 12 g. Including the heavier scalar mesons is only possible for smaller 
M. 

Finally, one should touch upon the convergence of the series (10) over n. For 
this aim we represent the coefficients ci(n) of (9) in the form (omitting the 
integrand) 



(lnZ)/27r 



c/(n) = c n (n) + c t2 (n) = 




M (lnZ)/2vr 



so that L(n) of (9) is equal to L\(n) + ^(n) respectively. Also it should be 
noted, as follows from the results of Ref . [4] , we have for barrier transparency 
T(k, I, n) of (9) that r(/x, l : n) — and at k — > +oo 

T(k,l,n) = \s u (k,l,n)\ 2 = 1 + 0(k~ v ) . (19) 



Further the asymptotic behaviour for large / and k « I 



\sn(k,l,n)\ = ^-(k-fi)[l + o(l)] (20) 

is crucial. To obtain it we should work using the equations (11)— (15) by con- 
ventional methods of mathematical analysis so that the detailed derivation of 
(20) lies somewhat out of the scope of the present paper and will be considered 
elesewhere. Using (20) we have the following chain of inequalities 



(ln0/2w | 12; +j; -21 {hll)/2lT 1 4, 

/ x f \s u rk + dk „ e Zl r , ,„ „ 9 ,ln I „ 
Ci(n)= j [HI — < B lir _/ A; 3 rfA; < B 2 e" 2 ^ < 5 3 e- 2 ' 



2Z 



with some constants = depending on /x so that (2Z + l)Qi(n) < -B/e 

with some constant B = B(/i) depending on /i. Since 



y Ze- 2, dZ = (|n| + l/2)e~ 2 ^/2 

= \n\ 
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then according to the Cauchy theorem 

oo 

Li(n) = A J2 (21 + l)Qi(n) < AB(\n\ + l/2) e - 2|n| . (21) 

l=\n\ 



At the same time due to (19) 

c l2 (n)<D J k + e~ 8nk+ dk ~ ^~JT ( 22 ) 



(lni)/27T 

with some constant D = D(/i) depending on ji. This entails 

(2i + l)lnZ j; ^lnH 

I" 



L 2 (n) ~ AD / ^ -/ dl ~ AD — . (23) 

7 r n 2 



Under this situation, according to (21) and (23), we obtain that L(n) ~ L 2 (n) 
and the series of (10) is convergent, i. e., L < +oo. Consequently, we can say 
that all the further computed luminosities are well defined and exist. 



3 Numerical results 

In view of (19), we have [denoting the Riemann zeta function as C( s )] 

J r{k,l,n)k+dk J kdk 1 1 

J e ^ k+ -l y e ^ k -l (8tt) 2U } 6-8 2 ' 1 } 

M V ' 

whilst the latter integral can be accurately evaluated using the trapezium 
formula on the [0, 3] interval. Having confined the range of k to the [/i, 3] 
interval, accordingly, it is actually enough to restrict oneself to < l,n <10- 
15 when computing L(n),L. We took < / < 14, < n < 10. Since the 
Wronskian of (11) does not depend on x, the latter should be chosen in the 
region where the potentials V + are already small enough. Besides, potentials 
V are really equal to when x <-12.5. We computed T(k,l,n) according to 
(11) at x — 500, where f ± were obtained from the Volterra integral equations, 
respectively, (12) and (16). For this aim, according to the methods developed 
for numerical solutions of integral equations (see e. g., Ref. [8]), those of (12) 
and (16) have been replaced by systems of linear algebraic equations which can 
be gained when calculating the conforming integrals by the trapezium formula, 
respectively, for the [-12.5, 500] and [500, 600] intervals. The sought values of 
were obtained as the solutions to the above linear systems. After this 
the derivatives (f^'j. were evaluated in accordance with (13) and (17) while 
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employing the values obtained for f^. The typical behaviour of T(k,l,n), for 
instance, for charged kaons is presented in Fig. 2. 

Fig. 2. Typical behavior of barrier transparencies for charged kaons. 

As an illustration Fig. 3 represents the untwisted L(0)/A and the all configura- 
tions L/A luminosities with A of (9) for charged kaons as functions of k while 
Fig. 4 presents luminosity for neutral pions (no interaction with monopoles). 
The areas under the curves give the corresponding values of L(0)/A and L/A. 

Fig. 3. Untwisted and all configurations luminosities for charged kaons. 

Fig. 4. Luminosity for neutral pions. 

In Table 1 the data on computation of the untwisted and the all configurations 
luminosities are represented for all meson species under consideration. Finally, 
we can introduce the total untwisted £ (as the sum over all species in the 
second column of Table 1) and the total £ (as the sum over all species in the 
third column of Table 1) luminosities multiplied by A that should be detected 
by an external observer near black hole. We shall have 

£ = 0.207299 • 10 23 erg • s" 1 , £ = 0.233520 • 10 23 erg • s" 1 



so that contribution owing to Dirac monopoles amounts to 11.2284 % of £. 
Table 1 

Untwisted and the all configurations luminosities 



Particle 


L(0)/A 


L/A 


Monopole contribution of L (%) 


7T° 


0.938559-KT 4 


0.938559-KT 4 





7T+ 


0.110489-KT 3 


0. 134436- 1(T 3 


17.8124 


7T~ 


0.110489-KT 3 


0.13443640" 3 


17.8124 


K° 


0.110660-KT 3 


0.110660-10- 3 





K° 


0.110660-KT 3 


0.110660T0- 3 





K+ 


0.110658-1CT 3 


0.13461640~ 3 


17.7979 


K~ 


0.110658-10- 3 


0.13461640- 3 


17.7979 
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4 Concluding remarks 



Although our algorithm works well only for small /i, it seems to us that includ- 
ing heavier mesons or considering more massive black holes will not essentially 
alter the results of the paper, so long as in either case /x will increase but it 
is clear that luminosity L — > at /i — > +00 and, as a consequence, the contri- 
butions from heavier mesons will be far smaller than the ones for pions and 
kaons. The more interesting case is that of spinor field which describes gen- 
uine elementary particles (electrons, neutrinos etc.) and, as has recently been 
shown [9], where also there exist TICs on black holes due to interaction with 
Dirac monopoles. We hope to shed some light on the spinor case elsewhere. 
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